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Abstract 

The novel dynamical features underlying soliton interactions in coupled nonlinear Schrodinger 
equations, which model multimode wave propagation under varied physical situations in nonlinear 
optics, are studied. In this paper, by explicitly constructing multisoliton solutions (upto four-soliton 
solutions) for two coupled and arbitrary iV-coupled nonlinear Schrodinger equations using the 
Hirota bilinearization method, we bring out clearly the various features underlying the fascinating 
shape changing (intensity redistribution) collisions of solitons, including changes in amplitudes, 
phases and relative separation distances, and the very many possibilities of energy redistributions 
among the modes of solitons. However in this multisoliton collision process the pair-wise collision 
nature is shown to be preserved in spite of the changes in the amplitudes and phases of the solitons. 
Detailed asymptotic analysis also shows that when solitons undergo multiple collisions, there exists 
the exciting possibility of shape restoration of atleast one soliton during interactions of more than 
two solitons represented by three and higher order soliton solutions. From application point of 
view, we have shown from the asymptotic expressions how the amplitude (intensity) redistribution 
can be written as a generalized linear fractional transformation for the ^-component case. Also we 
indicate how the multisolitons can be reinterpreted as various logic gates for suitable choices of the 
soliton parameters, leading to possible multistate logic. In addition, we point out that the various 
recently studied partially coherent solitons are just special cases of the bright soliton solutions 
exhibiting shape changing collisions, thereby explaining their variable profile and shape variation 
in collision process. 

PACS numbers: 42.81.Dp, 42.65.Tg, 05.45.Yv 
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I. INTRODUCTION 



The study of coupled nonlinear Schrodinger (CNLS) equations is receiving a great deal of 
attention in recent years due to their appearance as modelling equations in diverse areas of 
physics like nonlinear optics [1], including optical communications [2], bio-physics [3], mul- 
ticomponent Bose-Einstein condensates at zero temperature [4], etc. To be specific, soliton 
type pulse propagation in multimode fibers [1] and in fiber arrays [5] is governed by a set of 
iV-CNLS equations which is not integrable in general. However, it becomes integrable for 
specific choice of parameters [6,7]. On the other hand, the recent studies on the coherent 
[8] and incoherent [9] beam propagation in photorefractive media, which can exhibit high 
nonlinearity with extremely low optical power, necessitate intense study of CNLS equations 
both integrable and nonintegrable. The first experimental observation of the so-called par- 
tially incoherent solitons with the excitation of a light bulb in a photorefractive medium 
[10] has made this study even more interesting. In this context of beam propagation in a 
Kerr like photorefractive medium, the governing equations are a set of iV-CNLS equations 



where qj is the envelope in the jth mode, z and t represent the normalized distance along the 
fiber and the retarded time respectively, in the context of soliton propagation in multimode 
fibers. In the case of fiber arrays qj corresponds to the jth core. Here 2/x gives the strength 
of the nonlinearity. In the framework of N self trapped mutually incoherent wavepackets 
propagation in Kerr-like photorefractive media [11,12], qj is the jth component of the beam, 
z and t represent the coordinates along the direction of propagation and the transverse 
coordinate, respectively. The interesting property of the iV-CNLS equations of the form (1) 
is that they are integrable equations and possess soliton solutions. 

It is obvious from Eq. (1) that for iV = 1 it corresponds to the standard envelope soliton 
possessing integrable nonlinear Schrodinger equation, governing intense optical pulse propa- 
gation through single mode optical fiber [1,14]. For N = 2 case, it reduces to the celebrated 
Manakov model [13] describing intense electromagnetic pulse propagation in birefringent 
fiber. Manakov himself has made a detailed asymptotic analysis of the inverse scattering 




A? 




(1) 



p=l 
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problem associated with the system (1) for AT = 2 and identified changes in the polariza- 
tion vector [13]. However no explicit two soliton expression was given there. Very recently, 
Radhakrishnan, Lakshmanan and Hietarinta have obtained the bright one- and two-soliton 
solutions for this case [15] and have revealed certain novel shape changing (intensity re- 
distribution) collision properties. These Manakov solitons have been observed recently in 
AlGaAs planar waveguides [16] and precisely this kind of energy exchanging (shape chang- 
ing) collisions has been experimentally demonstrated in [17]. The results of Ref. [15] have 
led Jakubowski, Squier and Steiglitz [18] to express the energy redistributions as linear 
fractional transformations so as to construct logic gates. Later, Steiglitz [19] explicitly con- 
structed such logic gates including the universal NAND gate, based on the shape changing 
collision property, and hence pointed out the possibility of designing an all optical computer 
equivalent to a Turing machine, at least in a mathematical sense. However, results are scarce 
for N > 2 case of Eq. (1) though they are of considerable physical importance as mentioned 
earlier. 

The novel shape changing collision property exhibited by the 2-CNLS equations, which 
has not been observed in general in any other simpler (1 + 1) dimensional integrable system, 
requires a detailed analysis to identify the various possibilities and the underlying potential 
technological applications. In a very recent letter [20], the present authors have studied the 
multicomponent N-CNLS equations and shown that shape changing collisions occur here also 
with more possibilities of energy redistribution. It has also been briefly pointed out that 
the much discussed partially coherent solitons (PCSs) [11,12], which are of variable shape, 
namely 2-PCS, 3-PCS,..., Af-PCS, are special cases of the 2-soliton, 3-soliton,..., A^-soliton 
solutions of the 2-CNLS, 3-CNLS,..., A^-CNLS equations, respectively. The understanding of 
variable shapes [11,12] of these recently experimentally observed partially coherent solitons 
[21] in photorefractive medium and their interesting collision behaviour will be fascilitated 
by obtaining the higher order soliton solutions of the 2-CNLS and the A^-CNLS (N> 2) 
equations. 

In this paper, we wish to undertake a detailed analysis of the dynamical features associ- 
ated with soliton interactions in multicomponent A^-CNLS equations. There exists numerous 
interesting phenomena which one has to pay attention in order to realize the full potential- 
ities of these equations and the underlying novel soliton dynamics. Some of the important 
aspects include the following among others: 
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1. Explicit expressions for multisoliton solutions in multicomponent CNLS equations use- 
ful for analysis of interactions (as against formal expressions). 

2. Novel soliton interactions involving shape changing collisions. 

3. Dependence of shape changes and relative separation distances on amplitudes of the 
colliding solitons. 

4. Identification of different possibilities of energy redistributions among the different 
modes of the soliton during collision and obtaining generalized linear fractional trans- 
formations. 

5. State restoring properties in multisoliton solutions. 

6. Existence of partially coherent solitons, stationary and moving, as special cases of the 
above multisoliton solutions. 

7. Identification of multisoliton solutions as logic gates in multicomponent CNLS equa- 
tions. 

The present paper will be essentially devoted to the understanding of multisoliton inter- 
actions in N-CNLS equations, and its application in constructing logic gates and in iden- 
tifying partially coherent solitons as special cases of multisoliton solutions. In particular, 
in the present paper, we will deduce explicit expressions for multisoliton solutions (upto 
four-soliton solution), which can be easily generalized to arbitrary soliton case, 2-CNLS and 
then for arbitrary iV-CNLS equations. To start with, we will briefly consider the two-soliton 
solution to bring out the shape changing nature of soliton collisions, which can be quantified 
in terms of generalized linear fractional transformations (LFTs), and identify the changes 
in amplitudes, phases and relative separation distances among the solitons by carrying out 
appropriate asymptotic analysis. However the standard (shape preserving) elastic collisions 
can occur for specific choice of soliton parameters (initial conditions). More interestingly, we 
also point out that when more than two solitons collide successively, say three solitons, there 
exists the exciting possibility of restoration of the shape of one of the three solitons leaving 
the other two undergoing shape changes and we prove that the underlying soliton interac- 
tion is pair-wise. We give explicit conditions for the shape restoring property. Extending 
this analysis, one can easily check that in an M-soliton collision, it is possible to restore 
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the states of (M — 2) solitons after collisions. Such possibilities lead to the construction of 
optical logic gates of different types and generalized linear fractional transformations, as we 
will show in this manuscript. 

This paper is organized as follows. In Sec. II we briefly present the bilinearization 
procedure for the iV-CNLS equations. Explicit multisoliton solutions (upto four) of the 
2-CNLS equations are obtained in Sec. III. Then the generalization of these multisoliton 
solutions to iV-CNLS equations is given in Sec. IV. The two soliton collision properties of 
2-CNLS and their generalization to iV-CNLS equations are studied in Sec. V. In Sec. VI, 
we present a systematic procedure to identify the intensity redistribution among N modes 
in terms of a generalized linear fractional transformation which is the precursor to develop 
logic gates without interconnecting discrete components [18]. The interesting feature of 
the higher order soliton solutions, namely the pair-wise nature of collision of solitons, and 
the shape restoration property of the state of one soliton only in a three soliton collision 
process are presented in Sec. VII. In Sec. VIII we introduce the possibility of looking at 
the bright soliton solutions as logic gates, as an alternate point of view. Then in Sec. IX 
we demonstrate explicitly that for specific choices of the parameters of the bright soliton 
solutions various PCSs reported in the literature result. The collision properties of PCSs 
and the salient features of multisoliton complexes are also discussed. Sec. X is allotted for 
conclusion. Also in Appendix A we present the explicit form of the four-soliton solution. 



II. BILINEARIZATION 

The set of iV-CNLS equations (1) has been found to be completely integrable [6,7] and 
admits exact bright soliton solutions. Their explicit forms can be obtained by using the 
Hirota's bilinearization method [22] , which is straightforward. Any of the other soliton pro- 
ducing methodologies in principle is equally applicable, however, this paper is not concerned 
with the relative merits of the various methodologies. 

To start with, we make the bilinearizing transformation (which can be identified system- 
atically from the Laurent expansion [6]) 

qj = V> J = ^,...,N, (2) 
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to Eq. (1). This results in the following set of bilinear equations, 

(iD z + D*)gV\f = 0, j = l,2,...,iV, (3a) 

TV 

Dtf-f = 2^9 (n) 9 {n) \ (3b) 

n=l 

where * denotes the complex conjugate, g^'s are complex functions, while f(z,t) is a real 
function and the Hirota's bilinear operators D z and D t are defined by 

The above set of equations can be solved by introducing the following power series expansions 
for gr^'s and /: 

g U) = Xg U) + X * g U) + ^ 3 = XX ^ ^ (4a) 

/ = l + X 2 /2 + X 4 /4 + -, (4b) 

where x is the formal expansion parameter. The resulting set of equations, after collecting 
the terms with the same power in x, can be solved recursively to obtain the forms of g^'s 
and /. Though a formal closed form solution of the iV-soliton expression of Eq. (1) as a 
ratio of two (N x N) determinants can be given [23], it becomes necessary to deduce the 
explicit expressions (which is nontrivial) in order to understand the interaction properties 
at least for the lower order solitons. In the next section we will only present the minimum 
details. 



III. MULTISOLITON SOLUTIONS FOR N = 2 CASE 

As a prelude to understand the nature of soliton solutions for arbitrary iV-CNLS equa- 
tions, we first present the bright one- and two- soliton solutions of Eq. (1) with N=2 
(Manakov) case as given in Ref. [15] and then extend the analysis to obtain the explicit 
higher order soliton solutions. 

A. One-soliton solution 

After restricting the power series expansion (4) as 

9® = X9?\ 3 = 1,% / = l + X 2 / 2 , (5) 
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and solving the resulting set of linear partial differential equations recursively, one can write 
down the explicit one-soliton solution as 











l-l 





l + e m+m+R \ A2 J cosh (r} 1R + f) ' 



(6) 



where 771 = A;i(t + fc), A,- = Ql 1/2 , j = 1,2, and e H = \ kl+k *y — -■ Note 



that this one-soliton solution is characterized by three arbitrary complex parameters a{\ 
af\ and k\. Here the amplitude of the soliton in the first and second components (modes) 
are given by k\ R Ai and kiRA 2 , respectively, subject to the condition \Ai\ 2 + |A 2 | 2 = -, 
while the soliton velocity in both the modes is given by 2ku. Here km and ku represent 
the real and imaginary parts of the complex parameter k\. The quantity 



log 



(fei+fcf) 2 



2fcl B " 2fel 

denotes the position of the soliton. 



B. Two-soliton solution 

The two-soliton solution of the integrable 2-CNLS system has been obtained in Ref. [15] 
after terminating the power series (4) as 

9 U) = \H\ ! ■ VW- J = 1,2, (7a) 
/ = 1 + X 2 h + X%, (7b) 

and again solving the resultant linear partial differential equations recursively . Then the 
explicit form of the two-soliton solution can be written as 

Q,C?) e »ji _|_ a ^ e m _|_ e vi+vl+V2+Si 3 _|_ e vi+m+V2+ s 2 3 
Qj = - ^ , J = 1,2, (8a) 

where 

D — 1 _)_ e vi+vt+Ri _|_ e vi+V2+ s o _|_ e Vi+V2+s^ _|_ e v2+V2+R2 _|_ e m+vl+m+vi+R3 _ ^g^^ 

In Eqs. (8), we have defined 

„. - k (t + ik z) r s ° - Kl2 r Rl - — c R * - K22 

fci + Ki + k* h + k\ 

Slj = (h ~ h)^?^! - a^K-n) s 2j = ( fc 2 - «i)(«2 J)/t i2 - a?^) 
(h + fc*)(fc* + k 2 ) ' 6 ~ {k 2 + k* 2 ){k 1 + k* 2 ) 

eR ' A = 7Z I uJn! I t*L I r^ («u«22-«i2«2i), (8c) 



*|2 
2l 
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and 

E2 (n) (n)* 

n=l a i <*l o ,0^ 

The above most general bright two-soliton solution is characterized by six arbitrary complex 
parameters ki, k 2 , ctf \ and a 2 \ j = 1,2, and it corresponds to the collision of two bright 
solitons. Note that in Ref. [15], 5n, ^12, #21, and 5 22 are called as 5i, 5[, 5 2 , and 5^, 
respectively. The redefined quantities (%'s, i,j = 1, 2, are now used for notational simplicity. 



C. Three-soliton solution 

The two-soliton solution itself is very difficult to derive and complicated to analyse [15]. 
So obtaining the three-soliton solution is a more laborious and tedious task. However we 
have successfully obtained the explicit form of the bright three-soliton solution also. In order 
to obtain the three-soliton solution of Eq. (1) for the N = 2 case we terminate the power 
series (4a) and (4b) as 

g U) = x9? + x 3 gi j) + x 5 gi j \ (9a) 

/ = 1 + X 2 f2 + X 4 h + X 6 /e, j = 1, 2. (9b) 

Substitution of (9) into the bilinear Eqs. (3a) and (3b) yields a set of linear partial differential 
equations at various powers of x ■ The three-soliton solution consistent with these equations 
is 



Qj = 



D 

e m+V* 2 + r l3+$43 _|_ e V3+ri3+Vl+SB3 _|_ g»73 +^3 +»72 +*6j _|_ e Vl+m+V3+S7j _|_ e m+V2+V3+Ssi 

+ D 

e Vl+V2+V*j+ S 9j _|_ e Vl+Vl+V2+V2+ r l3+Tlj _|_ e Vl+V*+V3+V3+1l2+V2j 

+ D 

e m+V2+V3+ri3+rii+r3j 

+ Y) ' J = 1 ' 2 ' (10a) 
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X _l_ e m+v*+Ri _|_ e V2+V2+ R 2 _|_ e m+V3+R:i _|_ e *7i+T?|+<Si _|_ e vl+m+5t 

e m+V3+$20 _|_ g^+^i+^O _|_ g»?2+r?3+<530 _|_ g»?5+^3+<5| _|_ g»71 +r?? +r? 2 +r?J +i?4 
_|_ e fl+i;+'i3+l)3+-R5 _|_ e m+Vi+V3+V3+R6 _|_ g^l+^+^+^+TlO _|_ e Vl+Vt+Va+V2+ T w 
_^ e V2+ri2+m+V3+T20 e m+V2+vt+V3+r2 _|_ g»73+»7|+»7l+^+-r30 _|_ g*?3+r?3 +^+»72 +r| 
_)_g'?l+^+'?2+T?5+r?3+??3+-R7_ (10b) 





+ ikiz),i = 1, 


2,3, 


(fcl 


- /c 2 )(tt? ) /«2i - 


- a^K U ) 






f k 2 ) 


(fci 


- k 2 ){a^ ] K 22 - 


- ai ] ^i2) 




(fci + fc 2 *)(fc 2 - 


ffc 2 *) 


(h 


- k 3 )(a^ ] K 33 - 


- 0^13) 




(h + k* 3 )(k* 3 - 


f fcl) 


(k 2 


- k 3 )(a { 2 o) K 3l - 


- a? ) K2i) 




(kt + k 2 )(kt + k 3 ) 


(h 


- k 2 ){a[ j) K 23 - 


- « 2 j) ki 3 ) 



= <5 2j = (fci - k^ja^K^ - a^Ku) 
(h + kl)(kl + k 3 ) 

j4 . = (k 2 - k 3 )(a { 2 3) K 32 - aj j) K 22 ) 
(k 2 + k*)(k* 2 + k 3 ) 

j 6j = (k 2 - k 3 )(a { i ] k 33 - ai ] k 23 ) 
(k* 3 + k 2 )(k* 3 + k 3 ) " : 
(fci - k^ja^Kw - a^Ku) 
(fc 1 + fc 2 *)(fc 2 * + fc 3 ) 



(10c) 



c 2 + fc*; 



(fc 2 -fci)(fc 3 -fci)(fc3-fc 2 )(fc 2 *-fci) 
(fc x * + fc 1 )(fc 1 * + fc 2 )(fc 1 * + fc 3 )(fc 2 * + fci)(fc 2 * + fc 2 )(fc 2 * + fc 3 ) 

x a[ 3 \n 2 iK, 32 - k 22 k 31 ) + a 2 } (K,i 2 K, 3 i - k 32 ku) + Oi£\KuK 22 - k 12 k 2 i) 

(fc 2 -fci)(fc 3 -fci)(fc 3 -fc 2 )(fc 3 -fci) 
(fc* + fci) (fci* + fc 2 )(fc* + fc 3 )(fc 3 * + fci)(fc 3 * + fc 2 )(fc 3 * + fc 3 ) 

x af\n 33 K 2l - k 31 k 23 ) + a { i\n 3l K 13 - k u k 33 ) + a^\K 23 K U - ki 3 k 2 i)] , 

(fc 2 -fci)(fc 3 -fci)(fc 3 -fc 2 )(fc^-fc 2 *) 
(fc 2 * + fci)(fc 2 * + fc 2 )(fc 2 * + fc 3 )(fc 3 * + h)(k* + k 2 )(k* + fc 3 ) 

®i\k 22 k 33 - k 23 k 32 ) + a { 2 3 \n 13 K 32 - k 33 k 12 ) + ai\K l2 K 23 - K 22 /ti 3 )] , 

(lOd) 
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and 



e Rm = 



e R 5 = 



h _l_ h* ' 



m— 1,2, 3, e' 



<5io 



«42 



(fc 2 — fci)(A;| — AjJ) 


(A;? + fci) 


(fc* + fc 2 )(fci + k* 2 ){k* 2 + fc 2 ) 
(fc 3 -fci)(^-fci*) 


(fc* + fci) 


(A;* + fc3)(fc* + fc 1 )(fc* + fc3) 
(fc3-fc2)(fc 3 *-fc 2 *) 


(fc* + fc 2 ) 


(fc^ + fc 3 )(fc3+fc2)(fc 3 + fc3) 
(fc 2 -fcl)(fc3-fc?) 


(fc; + fci) 


(fc* + fc 2 )(fc 3 * + fci)(fc 3 * + fc 2 ) 

(fcl-fc 2 )(fc3-fc 2 ) 


(fc 2 + fci) 


(A£ + fc 2 )(fc^ + fci)(fc^ + fc 2 ) 
(fc 3 -fcl)(fc 3 *-fc 2 ) 


(fc* + fco 


(fc 2 * + fc 3 )(fc 3 +fcl)(fc 3 + fc 3 ) 



fci + fc^' 



[«33«11 
[/t 22 K 33 

[«11«23 
[K 22 Kl 3 



«13 
" fci + fc 3 

K l 2 K 2 l] 
«13«3l] 
K 23 K 32 ] 
«21«13] 
^1 2 ^ 2 3] 



K 23 



fc 2 + fc 3 *' 



I fci - fc 2 | 2 |fc 2 - fc3| 2 |fc3 



[«33«12 - «13«32] 
"fcll 2 



(fci + fc*)(fc 2 + fc 2 *)(fc 3 + fc 3 *)|fcx + fc 2 *| 2 |fc 2 + fc 3 *| 2 |fc 3 + fc*| 2 
X [(«11«22«33 - «11«23«32) + («12«23«31 - «12«21«33) 



+ («21«13«32 - «22«13«3l) 



E2 (n) (n)* 

n=l a i a i 



, i,Z = 1,2,3. 



(lOe) 



(lOf) 



(*« + fc;) 

The above three-soliton solution represents three soliton interaction in the 2-CNLS equations 
and is characterized by nine arbitrary complex parameters a-^'s and fcj's, i — 1, 2, 3, j — 1, 2. 
One can also check that the above general three-soliton solution of the 2-CNLS equations 
reduces to that of the solution given in Ref. [24] for the particular case = 1. Further, 
the form in which we have presented the solution eases the complexity in generalizing the 
solution to multicomponent case as well as to higher order soliton soutions. 



D. Four-soliton solution 



The expression is quite lengthy, but it is written explicitly in terms of exponential func- 
tions so as to check the pair-wise nature of collisions. We indicate the form in Appendix 
A. One can generalize these expressions for the arbitrary N case also. However, it is too 
complicated to present the explicit form and so we desist from doing so. 



11 



IV. MULTISOLITON SOLUTIONS FOR THE N-CNLS EQUATIONS 



As mentioned in the introduction, results are scarce for Eq. (1) with N > 2 and there 
exists a large class of physical systems in which the A-CNLS equations occur naturally. 
Further, in the context of spatial solitons in photorefractive media, each fundamental soliton 
can be "spread out" into several incoherent components [25], as defined by the polarization 
vectors. Obtaining one-, two-, and higher order soliton solutions of A-CNLS equations will 
be of considerable significance in these topics. In order to study the solution properties of 
such systems we consider the integrable A-CNLS equations (1). Following the procedure 
mentioned in the earlier section we obtain the one-, two-, and three- (as well as four-) soliton 
solutions of A-CNLS equations as given below. Particularly the so-called partially coherent 
solitons will turn out to be special cases of these soliton solutions ( see Sec. IX below). 



A. One-soliton solution 



The one-soliton solution of Eq. (1) is obtained as 

(gi, ga, ... , q N ) T = k 1R e irilI sech (v 1R + (A ± , A 2 , . . . , A N ) T , (11) 

where Vl = h(t + ik lZ ), A, = a? /A, A = (/i(E7=i l«Pl 2 )) 1/2 , e R = A 2 /(h + A;*) 2 , a® 
and fci, j=l, 2 . . . , N, are (N + 1) arbitrary complex parameters. Further kmAj gives the 
amplitude of the jih mode (j = 1,2,..., N) and 2ku is the soliton velocity in all the iV 
modes. 



B. Two-soliton solution 

The two-soliton solution of Eq. (1) can be obtained by following the procedure given for 
the 2-component case. It can be written as 

Qj = ~ ' ^ = 1 ' 2 '-'^ ( 12 ) 

where the denominator D and the co-efficients e Rl , e R2 , e Rs , e 5 ° , e 5 o , e 5lj , and e &2j , bear the 
same form as given in (8c) and (8d), except that j now runs from 1 to N and that k^'s are 
redefined as 

EN (n) (n)* 
- n=1 ^ 011 tl-12 (13) 
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One may also note that the above two-soliton solution depends on 2(N+1) arbitrary complex 
parmeters af \ k\, and k 2 , j = 1,2, N. 



C. Three-soliton solution 



Following the procedure given in the previous section we obtain the 3-soliton solution to 
the iV-CNLS equations as 

ti)„m i „,ti)„r)2 _|_ a ^ e V3 _|_ e m+Vi+m+Sij _|_ e vi+vt+V3+S2j _|_ e V2+V2+vi+s 3j 



a] e'' 1 + a, , e 

1j = 



D 

e m+V2+n.i+$43 _|_ e m+V3+m+^5j _|_ g^a+^+^+^j _|_ e v$+m+V3+S7j _|_ gm+^+^s+^j 



+ D 

+ ^ , J = 1,2,..., iV. (14a) 

Here also the denominator D and all the other quantities are the same as those given under 
Eq. (10) except for the redefinition of Kill S clS 

v-^TV (n) (n)* 

Ku = (k + kt) 1 ' ^ / = 1 ' 2 ' 3 - ( 14b ) 

It can be observed from the above expression that as the number of solitons increases the 
complexity also increases and the present three-soliton solution is characterized by 3(iV+ 1) 
complex parameters af \ ot% \ j = 1,2, ...,N, k\, k 2 and k^. 

The above procedure can be generalized to obtain four-soliton solution and higher order 
soliton solutions as discussed in the case of 2-CNLS equations straightforwardly and one 
can identify the iV-soliton solution of iV-CNLS will be dependent on N(N + 1) arbitrary 
complex parameters. 



V. SHAPE CHANGING NATURE OF SOLITON INTERACTIONS AND INTEN- 
SITY REDISTRIBUTIONS 

The remarkable fact about the above bright soliton solutions of the integrable CNLS 
system is that they exhibit fascinating shape changing (intensity redistribution / energy 
exchange) collisions as we will see below. This interesting behaviour has been reported in 
Ref. [15] for the two-soliton solution of the 2-CNLS equations. In a very recent letter [20], the 
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present authors have constructed the two-soliton solution of the 3-CNLS and generalized it to 
iV-CNLS, for arbitrary N, and briefly indicated similar shape changing collision dynamics of 
two interacting bright solitons. As these iV-CNLS equations arise in diverse areas of physics 
as mentioned in the introduction, it is of interest to analyse the interaction properties of 
the soliton solutions of 2-, 3-, and iV-CNLS equations. The collision dynamics can be well 
understood by making appropriate asymptotic analysis to the soliton solutions given in the 
previous sections. Such an analysis will then be used to identify suitable generalized linear 
fractional transformations in the next section, to obtain possible multistate logic. 



A. Asymptotic analysis of two-soliton solution of 2-CNLS equations 

To start with we shall briefly review the collision properties associated with the two- 
soliton solution (8) of the 2-CNLS equations discussed in Ref. [15] in order to extend the 
ideas to the iV-CNLS case. Without loss of generality, we assume that kj R >0 and ku>k 2I , 
kj = kj R + ikji, j = 1, 2, which corresponds to a head-on collision of the solitons (For the 
case ku = k 2 i, see Sec. IX below). For the above parametric choice, the variables r/jRS 
(real part of rjj) for the two solitons behave asymptotically as (i)r)w ~ 0, r] 2R — > ±oo as 
z — > ±oo and (n)i] 2R ~ 0, r]i R — > =poo as z — > ±oo. This leads to the following asymptotic 
forms for the two-soliton solution. (For other choices of km and ku, i — 1,2, similar analysis 
as given below can be performed straightforwardly), 
(i) Before Collision (limit z — > — oo) 
(a) Soliton 1 (r] 1R ss 0, r] 2R -^—oo): 



j fcifle^'sech (vir + y ) > ( 15a ) 




where 



1 2 

The quantity e Rl is defined in Eq. (8c). 



a 



(2) (h + kt) 



(15b) 



0, ViR ~ 


^oo): 


')- 




Q2 J 





k 2R e^'sech ( r, 2R + {R3 - Rl) ) , (16a) 
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where 

-(Ri+R 3 )/2 







( e 5 " \ 









(16b) 



The quantities in the above expression are again defined in Eq. (8c). 
(ii) After Collision (limit z — > oo) 

Similarly, for z — > oo, we have the following forms for solitons S\ and S 2 - 
(a) Soliton 1 (r]i R ss 0, i] 2R — > 00): 

Al+ \ 1R e^sech(r 1lR + iR3 - R2) ). (17a) 




where 



A\ + \ ( e &21 \ e ~(R2+R 3 )/2 
A\ + ) = \ e s ™ J (** + k l) ' 

(b) Soliton 2 (i] 2R ~ 0, 77^ ^—00): 





hi 


< A? \ 









where 



A\+ j \ «?> / (fe + K) 



(17b) 



W^sech U 2ii + ^), (18a) 



(18b) 



In the above expressions for Si and 5*2 after collision the quantities e R2 , e Rz , e S21 and e 522 
are defined in Eq. (8c). 



B. Asymptotic analysis of the 2-soliton solution of N-CNLS equations 

We require the asymptotic forms of the 2-soliton solutions for arbitrary N case in the 
following section in order to identify a generalized linear fractional transformation for the 
amplitude redistribution among the components. To get the asymptotic forms of 2-soliton 
solution of the A-CNLS case, as may be checked by a careful asymptotic analysis along the 
lines of the N = 2 case, we simply increase the number of components in the A ± vectors 
above up to N (A ± = (Af, Af, ^) T ) by adding two more complex parameters atf\ a 2 \ 
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i = 3,4, ...,N, to each of the components so that the forms of the quantities e Rl , e R2 , e Ri , 
e Su ^ e <5i2 ? e S2i ^ e &22 j n (g c ) rema i n the same as above except for the replacement of the 
range of the summation in n a (Eq. (8d)) from n = 1, 2 to n = 1, 2, N. As an example, in 
the following we give the asymptotic forms of two-soliton solution of the iV-CNLS equations 
with N = 3, for the case k [R > 0, I — 1, 2, and k u > k 2 i- For other possibilities similar 
analysis can be made, 
(i) Before Collision (limit z — > — oo) 
(a) Soliton 1 (r)i R ~ 0, ^r^—oo): 



(a\-\ 

A\- 



k 1R e lVu sech r} 1R + 



Ri 



(19a) 



where 



a 



(2) 



(3) 



(fcl + fcj)' 



(19b) 



(b) Soliton 2 (i] 2 r w 0, 77^ 



52 



oo 

/ 



2-\ 



k 2R e iV2I sech r] 2R + 



(R 3 - Ri) 



(20a) 



where 



/ e «5ii \ 



5l3 



e -(fll+fl3)/2 

(fe + A; 2 *) • 



(20b) 



(ii) After Collision (limit z — > oo) 

(a) Soliton 1 (rji R ~ 0, r/ 2 i? — > oo) 



52 

V 5 3 y 



/ci R e ir?1/ sech 77^ + 



(i?3 - #2) 



(21a) 
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where 



A? 



(b) Soliton 2 (i] 2 r w 0, ^ ->- oo): 



\ efe V 



,-(fl 2 +fl3)/2 



5i 

52 

V 5 3 y 



4 2 + 



where 



( „<U \ 



a 2 

V«2 / 



-fl 2 /2 



(fc 2 + fc 2 *) 



(21b) 



(22a) 



In the above expressions, the forms of the quantities e Rj , e Sij , i = 1 , 2 , j 
identified from Eqs. (12) and (13) with N = 3. 



(22b) 



1, 2, 3, can be 



i. Intensity redistribution 

The above analysis clearly shows that due to the interaction between two copropagating 
solitons S 1 and S 2 in an iV-CNLS system, their amplitudes change from A^~k 1R and A 2 ~k 2 R 
to A^kiR and A 2+ k 2R , j — 1, 2, . . . , TV, respectively. However, during the interaction process 
the total energy of each of the solitons is conserved, that is 

f>ri 2 =x>rM- (23) 

3=1 3=1 ^ 

Note that this is a consequence of the conservation of L 2 norm. Another noticeable observa- 
tion of this interaction process is that one can observe from the equation of motion(l) itself, 
that the intensity of each of the modes is separately conserved, that is, 

/CO 
\qj\ 2 dz = constant, j = 1, 2, . . . , N. (24) 
-oo 

The above two equations (23) and (24) ensure that in a two soliton collision process (as well 
as in multisoliton collision processes as will be seen later on) the total intensity of individual 
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solitons in all the N modes are conserved along with conservation of intensity of individual 
modes (even while allowing an intensity redistribution). This is a striking feature of the 
integrable nature of the multicomponent CNLS equations (1). The change in the amplitude 
of each of the solitons in the jth mode can be obtained by introducing the transition matrix 



TL j = 1, 2, N, 1 = 1, 2, such that 



A? 



(25a) 



The form of T-'s can be obtained from the above asymptotic analysis as 



1 [ a*Jy « 12 



1 - A 



2 1 



VI - AiA 2 



, j = l,2,...,N, 



(25b) 



where 



a,2 = (k 2 + hi) 



(n) (n)* 
™2 



and 



a i 



«21 



N 

(ki — k 2 ) a^a 

n=l 



yl — A1A2 



1/2 



(25c) 



1-Ai 



in which 



ai = (h + kl) 



N 



1,2, 



(25d) 



(ki - k 2 ) ^2 a i 



(n)* (n) 



Q.-, 



n=l 



1/2 



(25e) 



In the above expressions Ai = ^ and A2 = where k^s , i,l = 1,2, are defined in Eq. 
(13). Then the intensity exchange in solitons Si and S 2 due to collision can be obtained by 
taking the absolute square of Eq. (25b) and (25d), respectively. 

The above expressions for the components of the transition matrix implies that in general 
there is a redistribution of the intensities in the iV modes of both the solitons after collision. 
Only for the special case 



a 



(2) 



a 



(N) 



a 



(2) 



a 



(26) 



there occurs the standard elastic collision. For all other choices of the parameters, shape 
changing (intensity redistribution) collision occurs. 
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The two conservation relations (23) and (24) allow the intensity redistribution to take 
place in definite ways. In general, for N-CNLS equations the intensity redistribution in a 
two soliton collision can occur in 2^ — 2 ways. Denoting E and S as enhancement and 
suppression respectively, either complete or partial, of the intensity of corresponding modes, 
we table below the possibilities of intensity redistribution for the case N = 2 and N — 3. 



Case 


Qi 


12 


1 


E 


S 


2 


S 


E 



(a) iV = 2 case 



Case 


Qi 


?2 


03 


1 


E 


S 


s 


2 


S 


E 


s 


3 


S 


S 


E 


4 


s 


E 


E 


5 


E 


S 


E 


6 


E 


E 


S 



(b) N = 3 case 

Table-I. Possible combinations of intensity redistribution among the modes of soliton Si in 

the two soliton collision process. 

For each of the above choices of Si, the form of S2 is determined by the conserved quantity 
(24) for the intensities of the individual modes. For illustrative purpose, we have shown in 
Fig. 1 and Fig. 2 a few of such possibilities of intensity switching for the N = 2 and N = 3 
cases, respectively. 

FIG. 1: Two distinct possibilities of the shape changing two soliton collision in the integrable 
2-CNLS system. The parameters are chosen as (a) k\ = 1 + i, ki = 2 — i, = = = 1, 
a W = 39i|0i ; (b) jfex = 1 + i, jfe 2 = 2 - i, = 0.02 + O.li, af ] = = = 1. 
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FIG. 2: Intensity profiles of the three modes of the two-soliton solution in a waveguide described 
by the 3-CNLS (Eq. (1) with N = 3) showing different dramatic scenarios of the shape changing 
collision for various choices of parameters. 



2. Phase shifts 



Further, from the asymptotic forms of the solitons Si and S 2 , it can be observed that 
the phases of solitons Si and S 2 also change during a collision process and that the phase 
shifts are now not only functions of the parameters ki and k 2 but also dependent on ck^'s, 
i — 1,2, j — 1,2, N. The phase shift suffered by the soliton Si during collision is 



(i? 3 — Ri — R 2 ) 



log 



\ki — k 2 \ 2 (KnK 22 — Ki 2 K 2 i) 

\ki + k 2 \ 2 KiiK 22 



where /%'s are defined in Eq. (13). Similarly the soliton S 2 suffers a phase shift 

(R 3 — R 2 — Ri) 



(27) 



(28) 



Then the absolute value of phase shift suffered by the two solitons is 



l$l = I* 1 ! 



(29) 



Let us consider the case N = 2. For a better understanding let us consider the pure 
elastic collision case (a^ : = otf* ■ offi) corresponding to parallel modes. Here the 
absolute phase shift (see Eq. (29)) can be obtained as 



l$l = 



log 



\h-k 2 \ 2 

\ki + k* 2 \ 2 



log 



\ki — k 2 



\ki + k* 2 \ 

Similarly for the case corresponding to orthogonal modes (a^ 
the absolute phase shift is found from Eqs. (27)- (29) to be 

\ki - k 2 \ 



4 1 ' 



(2) 

oo, a \ 



(2) 



m = 



log 



Ifcl 



^2 1 



(30) 
= 0) 

(31) 



The absolute value of the phase shift takes intermediate values for other choices of the 



parameters a-^'s, i — 1,2, j — 1,2, ...,7V. Thus phase shifts do vary depending on a\ J>, s 
(amplitudes) for fixed fcj's. In Fig. 3, we plot the change of |<E»| as a function of when 
it is real, at ct[^ = a 2 ^ = 1, = 39 g 9 80t , k± — 1 + i, and k 2 = 2 — i. Similar analysis can 
be done for the N = 3 case and for the arbitrary iV case. 



Cj)> 
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FIG. 3: Plot of the magnitude of phase shift as a function of the parameter , when it is real (for 
illustrative purpose), see Eqs. (29-31). The other parameters are chosen as k\ = 1 + i, &2 = 2 — i, 
a f) = 4 2 ) = 1 and 4 ] = ™. 

3. Relative separation distance 

Ultimately the above phase shifts make the relative separation distance tf 2 between the 
solitons (that is, the position of S 2 (at z — > ±oo) minus position of Si (at z — > ±oo)) also 
to vary during collision, depending upon the amplitudes of the modes. The change in the 
relative separation distance is found to be 



+ _ 

K1RK2R 

Thus as a whole the intensity profiles of the two solitons in different modes as well as the 
phases and hence the relative separat: 
vary as a result of soliton interaction. 



At 12 = r 12 - t+ = ^\ ( 32 ) 



phases and hence the relative separation distance are nontrivially dependent on a^'s and 



VI. GENERALIZED LINEAR FRACTIONAL TRANSFORMATIONS AND MUL- 
TISTATE LOGIC 

The intensity redistribution was characterized by the transition matrix as given in Eq. 
(25) in the previous section. Interestingly, this redistribution can also be viewed as a linear 
fractional transformation (LFT) as already pointed out by Jakubowski et al. [18]. However, 
no systematic derivation of such a connection was made. In this section, we point out 
that in fact a reformulation of Eq. (25) allows one to deduce such a LFT in a systematic 
way. This in turn allows us to generalize the procedure to the iV-component case leading 
to a generalized LFT for the amplitude change during soliton collision thereby leading to a 
multistate logic. 



A. N=2 case 



For the N = 2 case, the amplitude change in the two modes of soliton 1 after interaction 
given by Eq. (25) can be reexpressed by the following transformation, which can be deduced 
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from comparison of expressions (15b) and (17b): 



A\ + = YC 2l A\- + YC 22 A\-. 



(33a) 



Here 



r — t(a\ , a\ , a\ , a\ 



(2)„,(2)>M 



«12 



^11^22 



-1/2 



(33b) 



in which a 2 is given in Eq. (25c). The forms of C^'s, i, j = 1, 2, read as 

C 12 = -^ 1) 4 2) *(A; 2 + ^), 

C 21 = -4 2) ^ 1} *(A; 2 + A; 2 *), 

C22 = af'a^^ + ^j + afa?^!-^). 



(33c) 
(33d) 
(33e) 
(33f) 



Note that the coefficients C^-'s are independent of a^'s and so of A[~ and A 2 ~, that is the 
a-parameters of soliton 1. Then from Eqs. (33a) the ratios of the A^'s, i,j = 1,2, can be 
connected through a linear fractional transformation(LFT). For example, for soliton 1, from 

(33a), 

(34) 



A\ + C llP \- 2 + C 



12 



C 2 ip\ >2 + C 22 



1- A 1 ^ 

where Pi j2 =^t^, in which the superscripts represent the underlying soliton and the subscripts 
represent the corresponding modes. The quantities p}t,, p\~ 2l Cn, C 12 , C 21 , C 22 , in Eq. (34) 

Pl) , h**£, h* and ((1 - h*)p L + j^j , respectively, 



are same as the quantities pn, pi, 



l-h* 

P*L 



given by Eq. (9) in Ref. [18] in an adhoc way. Thus the state of S\ before and after 
interaction is characterized by p\~ 2 and p}t,, respectively. It is to be noticed that during 
collision fcj's, % — 1, 2, are unaltered. The LFT has been profitably used in Ref. [19] to 
construct logic gates, associated with the binary logic p = [0,1]. Similar analysis can be 
done for the soliton 2 also. 
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B. N=3 case 



Extending the above analysis, straightforwardly one can relate the A^'s, j = 1, 2, 3, for 



soliton 1, from Eqs. (19b) and (21b), as 

A\ + = VC^A 1 - + YC l2 A\- + TC 13 A l f, 

A\ + = YC 2l A\- + YC 22 A\- + TC 23 Al~, 

A\ + = rC 31 A{- + YC 32 A\~ + TC 33 A\-, 

where 



(35a) 
(35b) 
(35c) 



r — t(a\ ,a\ , a\ , a\ ,a 2 , a 



(2)^,(2)* 



X 



1/2 



(2)^,(2)* 



„( 3 )«( 3 )* 



1^12 1 ^11^22. 

in which g^'s are redefined as 



a 2 = (jfe 2 + kl)[(ki - k 2 )(a?aP* + af ] a { 2 2) * + a^af*)] 1 ' 2 



(35d) 



(35e) 



and Kus can be written from Eq. (13) with N = 3. Note that the form of T is a straightfor- 
ward extension of the N = 2 case. In the above equations the coefficients C^-'s, i, j = 1, 2, 3, 
for the 3-CNLS case can be written down straightforwardly by generalizing the expressions 
(33) corresponding to the two-soliton solution of the two component case. 

Thus in the two soliton collision process of the N = 3 case, for soliton 1 we obtain the 
generalized Mobius transformation, 

i+ _ A\+ _ Cnpira + C^ + Ci 



-13 



A\ + C 3 ip\ 3 + C 32 p\ 3 + c 33 
i+ A^_ _ C 21 p{~ 3 + C 22 p\~ 3 + C 23 

A 3 + C3lPl,3 + C*32P2,3 + C33 



(36a) 
(36b) 



where p{ 3 = and p\ 



,3 



Similar relations can be obtained for the soliton 2 also. 



C. Arbitrary N case 

Proceeding in a similar fashion one can construct for the soliton Si a generalized linear 
fractional transformation for the N-component case also which relates the p vectors before 
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and after collision. 




(37a) 



A)f J2f=i c Njp] :N 



with the condition 



Pnn 



= 1. 



(37b) 



Here p\ N 



Jr. Similar expression can be obtained for soliton 2 also. 



The above generalization paves way not only for writing down the bilinear transformation 
but also to identify multistate logic. For example, in the iV = 3 case, the following states are 
possible: p = [pi, P2] = [(0, 0), (0, 1), (1, 0), (1, 1)], where the logical '0' state can stand for 
the complex valued p state corresponding to a suppression of the amplitude in that mode, 
while the logical '1' state may correspond to enhancement (including no change), which can 
be used to perform logical operations, whereas in the N = 2 case we have only the two state 
logic, p = [0, 1]. This shows that for N > 2, we will get multistate logic and we believe that 
such states can be of distinct advantage in computation. This kind of study is in progress. 

VII. HIGHER ORDER SOLITON SOLUTIONS AND THEIR INTERACTIONS 

Now it is of interest to study the nature of multisoliton collisions making use of the 
explicit forms of multisoliton solutions given in Sees. Ill and IV. Due to the complicated 
nature of the above bright soliton expressions, it becomes nontrivial to identify the nature 
of the collision process. In his paper [13], Manakov pointed out that in general an N-soliton 
collision does not reduce to a pair collision due to the nontrivial dependence of the amplitude 
of a particular soliton before interaction on the other soliton parameters. In this section by a 
careful asymptotic analysis of the 3-soliton solution (10) of the 2-CNLS equations, which can 
be deduced to the A-CNLS case without any difficulty, we explicitly demonstrate that the 
collision process indeed can be considered to occur pair-wise and thereby making Manakov's 
statement in proper perspective and clearer. One can carry out a similar analysis for the 
four-soliton solution given in Appendix A, generalizing which one can show that in the higher 
order solitons of CNLS equations also the collision is pair-wise. Such an analysis also reveals 
the many possibilities for energy exchange among the modes of the solitons, including the 
exciting possibility of state restoration in higher order soliton solutions, a precursor to the 
construction of logic gates. 
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A. Asymptotic analysis of 3-soliton solution of 2-CNLS equations 



Considering the explicit three soliton expression (10), without loss of generality, we as- 
sume that the quantities k 1R , k 2R , and k 3R are positive and k\j > k 2 j > k 3I (For the equal 
sign cases ku = k 2 i = k 3 j, see Sec. IX below). One can carry out similar analysis for other 
possibilities of kn's , i = 1,2,3, also as discussed below. Then for the above condition the 
variables ^_r's, % = 1,2,3, for the three solitons (Si, S 2 , and S 3 ) take the following values 
asymptotically: 

(i) r} 1R ~ 0, T] 2R -> ±oo, i] 3R -> ±oo, as z ^ ±oo, 

(ii) r] 2R w 0, 7] 1R -> =foo, 7] 3R ->■ ±oo, asz^ ±oo, 

(hi) T] 3R W 0, Tj 1R -»■ =FOO, 7? 2 R -> ±oo, as2^ ±oo. 

Defining the various quantities R^s, i = 1, 2, 7 , / = 1, 2, 9, j = 1,2, r mj -'s, and 
r m0 's, m = 1,2,3, as in Eq. (10) we have the following limiting forms of the three-soliton 
solution Eq. (10). 
(i)Before Collision (limit z — > — oo) 



(a) Soliton 1 (r) 1R ps 0, ^ — > -oo, i] 3R -> -oo): 



*> 


H 


















l-l 


U 2) 



k 1R sech [r) 1R + ^)e 



e 2 



(b) Soliton 2 (rj 2R ps 0, 



(c) Soliton 3 (rj 3R ps 0, 77^ -> oo, ?/ 2/? 






H 




92 J 








l-l 


f e T " 






1 e r,J 



(fci + fef)' 
— oo): 

/c 2j RSech ^T}2R + 

(R\+Rj) 

e 2 
oo): 

A^sech ^ 3i? ± 

(fl 4 +fl 7 ) 
e 2 

(*3 + *5) ' 
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(38a) 
(38b) 



i?4 — R\ 



(39a) 
(39b) 



(40a) 
(40b) 



(ii)After Collision (limit z — > +00) 

(a) Soliton 1 (i] 1R m 0, i] 2 r -> 00, i] 3R -> 00): 




Tin * * - 



i? 7 — R 6 




)-l 




(«6+«7) 
I" 2 • 




ie I!2 J 


(h + kt) 


? -> -00, ?73 R 


— > 00): 


)-l 




| £; 2jR sech ^ 2jR + 2 


)-l 


( e Stn \ 


(H 3 +flfi) 




> e <5 62 J 





(c) 5ofefora 5 (ry 3 R ~ 0, 77^ -> -00, r? 2i? -> -00): 





1-1 


















1-1 





/c 3/? sech ( ?73_r + — < < 



"73/ 



_^3 

e 2 



w (fc 3 + fcs) 



(41a) 
(41b) 



(42a) 
(42b) 



(43a) 
(43b) 



B. Transition elements 

The above analysis clearly shows that during the three soliton interaction process, there 
is a redistribution of intensities among these solitons in the two modes along with amplitude 
dependent phase shifts as in the case of the two soliton interaction. The amplitude changes 



can be expressed in terms of a transition matrix Tj as 

4+ = Tj4-,j = l,2, 1 = 1,2,3. 



(44) 



Explicit forms of the entries of the transition matrix quantifying the amount of intensity 
redistribution for the three solitons are as follows. 
Soliton 1: 




(45a) 
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Soliton 2: 



T 2 J I e «5 62 -5 12 ' 



Soliton 3: 



T l\_l a 3 )e ^ \ -(R3-R4-R7) 



e ' 2 . (45c) 
The various quantities found in the above equations are defined in Eq. (10). 

C. Phase shifts 

Now let us look into the phase shifts suffered by each of the solitons during collision. 
These can be written as 

5 1 : S 1 = R i- R 2 *- R \ (46a) 

5 2 : $ 2 = R s- R *- R ± + R \ (46b) 

5 3 : $ 3 = ^-^ + it>4 . (46c) 

Here the quantities R\, R2,..., R7 are as given in Eq. (10). Note that each of the phase 
shifts $ 1 , $ 2 , and $ 3 contains a part which depends purely on fcj's, i = 1,2,3, and another 
part which depends on the amplitude (polarization) parameters a-^'s along with k^s. 

D. Relative separation distances 

As a consequence of the above amplitude dependent phase shifts, the relative separation 
distances between the solitons tfj (position of Sj (at z — > ± 00) - position of Si (at z — > ± 
00), i ^ j, i < j, i,j = 1, 2, 3) also varies as a function of amplitude parameters. The change 
in the relative separation distances (Atij = t~j — t^) can be obtained from the asymptotic 
expressions (38-43). They are found to be 

&k2R - $ 2 k 1R 
Ati2 = : : , (47a) 

Atis = ( 47b ) 

kiRk^R 

At 2 3 = ; — ; , (47c) 

kiRk-zR 
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where $ J 's, j = 1, 2, 3, are defined in Eq. (46) and kj^s represent the real parts of k/s. 



E. Nature of Collision 

Now it is of interest to look into the nature of the collisions in the three soliton interaction 
process, that is, whether it is pair-wise or not. This can be answered from the asymptotic 
expressions presented in Eqs. (38) to (46). For example, let us consider soliton 1 (Si). 
The net change in the amplitudes of the two modes of soliton Si is given by the transition 
amplitudes T/, % = 1, 2, that is, 



i+ 




(48) 



where T± and T\ are defined in Eq. (45a). The above form of transition relations is obtained 
by expanding (44). 

Let us presume first that the collision process is a pair-wise one and then verify this 
assertion. According to our assumption ku > k 2 i > k 3I , and so the first collision occurs 
between Si and S 2 as shown schematically in Fig. 4. Then during collision with S 2 , the two 

FIG. 4: A schematic three soliton collision process (for the choice km, k 2 R, k^R > 0, ku > k 2 j > 
%r). The effect of phase shifts are not included in the figure. 

modes of Si change their amplitudes (intensities) by T\ and T 2) respectively. Their forms 
follow from the basic two soliton interaction process discussed in Sec. V, Eq. (25b). This 
can be expressed in mathematical form as 



where 



ft 




(49a) 



rpl J \ e°32 

J 2 / \ ~m 



e"3i 

(i) I -(fl 4 +fl 2 -fii) 

" ' " " (49b) 



2 



a. 



Again the above expressions can be obtained straightforwardly from Eq. (25b) with N = 2. 

Now the resulting soliton (Si), after the first collision, is allowed to collide with the third 
soliton (S 3 ) (see Fig. 4). From the asymptotic expressions (38-45) and using the above Eqs. 
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(49), it can be shown that 



1+ 




(50a) 



where 



3 T31— 531 



-(R G + R 7 -R 4 -R 2 ) 
2 



f\ I \ e T ' i2 ~^ i2 



(50b) 



However using Eq. (49) in (50a), we can write 




(51a) 



(51b) 

If this is the collision scenario, then the right hand sides of Eq. (48) and Eq. (51b) should 
be the same, that is, 

Tl = flfl (51c) 
Tl = f\f\. (51d) 

This can be easily verified to be true directly from the expressions (45) and (49-50). In 
a similar fashion, for the other two solitons also the transition matrix can be shown as a 
product of two matrices corresponding to two collisions, respectively. 

Now let us look at the phase shifts. It is also of necessary to identify whether the total 
phase shift acquired by each soliton during the three soliton collision process is a result 
of two consecutive pair-wise collisions or not. In this regard, we again focus our attention 
on soliton 1 (Si) first. Let us assume the collision to be pair-wise. Then one can write 
the phase shift suffered by Si during the collision based on the analysis of the two soliton 
collision process. Following Eq. (27) (with appropriately changed notations) we can write 
the expression for the phase shift suffered by S\ on its collision with S% as 

6 = R '~ R ^- R \ (52) 

Now the outcoming form of Si (which is Si) is allowed to interact with S3 (see Fig. 4). The 
phase shift during this second collision can again be found from the asymptotic expressions 
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(38-43) as 



- R 7 — R 6 — R4 + R 2 , s 

d = 2 • ( 53 ^ 



On the other hand, from the asymptotic expressions (46), the total phase shift suffered by 
Si in a three soliton collision process can be written as 

5 = R ?- R ^- R i (54) 
= 5 + 5. (55) 

Thus the total phase shift suffered by the soliton 1 is the sum of the phase shifts suffered by 
it during pair-wise collisions with soliton 2 and soliton 3, respectively. Similar conclusions 
can also be drawn on the phase shifts suffered by the other two solitons as well. Thus the 
above analysis on the changes in the amplitudes and phase shifts during the three soliton 
collision process establishes the fact that the collisions indeed occur pair-wise. 

It may be noted that the above results also imply that the three soliton collision process 
is associative and independent of the sequence in which collisions occur, that is whether 
the collision occurs in the order £1 — > S 2 — > S 3 or Si — > S 5 — > S 2 . This property has been 
anticipated in the numerical study of Lewis et al. [26] , which is now rigorously proved here. 



F. Intensity redistributions and shape restoration 

The asymptotic analysis not only explains the nature of the collision process, but also 
characterizes the collision process. It is clear from the above analysis of the three-soliton 
solution that in general there is an intensity redistribution among the three solitons due to 
pair-wise interaction in all the two modes along with amplitude dependent phase shifts as 
in the two soliton interaction, subject to conservation laws. We have analysed the various 
three soliton collision scenarios below. 



1. Elastic collision 

The standard elastic collision property of solitons results for the special case : ■ 
= oif* : ■ otf\ The magnitude of the transition elements |T-|, j = 1,2, and I = 
1, 2, 3, becomes one for this choice of parameters and there occurs no intensity redistribution 
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FIG. 5: Intensity profiles \qi\ 2 and \q2\ 2 of the two modes of the three-soliton solution of the 2-CNLS 
equations, representing elastic collision, with the parameters chosen as k\ = 1 + i, k 2 = 1.5 — 0.5?, 
k 3 = 2 - i, a? = a« = a« = af ] = = a< 2) = 1. 

FIG. 6: Intensity profiles \q\\ 2 and \q 2 \ 2 of the two modes of the three-soliton solution of the 2- 
CNLS equations, representing the shape changing (intensity redistribution) collision process for 
the choice of the parameters, k± = 1 + i, k 2 = 1.5 - 0.5i, k 3 = 2 — i, a>P = 39 ~ 9 80 ' , (dp = 39 g 9 8(h , 
a { p = 0.3 + 0.2s, a { p = 0.39, a { p = a { p = 1. 

among the modes except for phase shifts. This is shown in Fig. 5 for the parametric choice 
a\ j) = 1, I = 1, 2, 3, j = 1, 2, jfex = 1 + i, k 2 = 1.5 - 0.5i, and A; 3 = 2.0 - i. 

Shape changing (intensity redistribution) collision 

For general values of the parameters ck^'s, there occurs shape changing collisions among 
the three solitons, however, leaving the total intensity of each of the solitons conserved, that 
is, l-A^I 2 + l-A^p = -, / = 1, 2, 3. This intensity redistribution is accompanied by amplitude 
dependent phase shifts and changes in the relative separation distances of the solitons as 
discussed above. They can be calculated from the expressions (38-47). One such shape 
changing interaction is depicted in Fig. 6 for illustrative purpose. The parameters chosen 
are fci = 1 + i, k 2 = 1.5 - 0.5i, k 3 = 2 - i, op = ^f^i, a { p = a { p = 0.3 + 0.2i, 

(2) (2) (2) 

a\ = 0.39, «2 — a 3 — 1 ■ In this figure we have shown the scenario in which the three 
solitons in the two modes have different amplitudes (intensities) after interaction when 
compared to the case before interaction. Here Si is allowed to interact with 5*2 first and 
then with S3. Due to this collision, in q 1 mode the intensity of Si is suppressed while that 
of S 2 is enhanced along with suppression of intensity in S 3 . On the other hand, the reverse 
scenario occurs in the q 2 mode for the three solitons Si, S 2 , and S3. 

3. Shape restoration of any one of the three solitons 

The asymptotic analysis also shows that there is a possibility for any one of the three 
solitons to restore its shape (amplitude/intensity) during collision. In this connection, let us 
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look into how the shape restoring property of Si occurs during its collision with the other 
two solitons (say S 2 and S3). We have already shown that the collision process is a pair-wise 
one. Then the three soliton collision process is equivalent to two pair-wise collisions. Let 
the first collision be parametrized by the parameters a ± \ a x \ cttp, ct^\ fci, and k 2 . Now 
we exploit the arbitrariness involved in choosing the parameters and ajp in the second 
collision process in order to make the net transition amplitude of Si to be unity, leaving the 
other two transition amplitudes of S 2 and S3 to vary, that is, 



This condition will make the soliton Si only to be unaffected at the end of the three soliton 
collision process. Then the equations corresponding to this condition are 



Air + A 2R x - A 2I y + A 3R (x 2 - y 2 ) - 2A 3I xy + A 4R x + A 4I y + A 5R (x 2 + y 2 ) 

+A m (x 3 + xy 2 ) - A eI (x 2 y + y 3 ) + A 7R (x 2 - y 2 ) 
+2A 7I xy + A 8R (x 3 + xy 2 ) + A 8I (x 2 y + y 3 ) + A 9R (x 2 + y 2 ) 2 = 0, (57a) 
An + A 2R y + A 2I x + 2A 3R xy + A 3I (x 2 - y 2 ) + A 4I x - A 4R y + A 5I (x 2 + y 2 ) 

+A 6I (x 3 + xy 2 ) + A 6R (x 2 y + y 3 ) - 2A 7R xy 
+A 7I (x 2 - y 2 ) - A 8R (x 2 y + y 3 ) + A 8I (x 3 + xy 2 ) + A 9I (x 2 + y 2 ) 2 = 0, (57b) 
Bin + B 2R x - B 2I y + B 3R (x 2 - y 2 ) - 2B 3I xy + B 4R x + B 4I y + B 5R (x 2 + y 2 ) 

+B 6R {x 3 + xy 2 ) - B 6I {x 2 y + y 3 ) + B 7R {x 2 - y 2 ) 
+2B 7I xy + B 8R (x 3 + xy 2 ) + B 8I {x 2 y + y 3 ) + B 9R (x 2 + y 2 ) 2 = 0, (57c) 
B u + B 2R y + B 2I x + 2B 3R xy + B 3I {x 2 - y 2 ) + B 4I x - B 4R y + B 5I (x 2 + y 2 ) 

+B 6I {x 3 + xy 2 ) + B m {x 2 y + y 3 ) - 2B 7R xy 
+B 7I (x 2 - y 2 ) - B 8R (x 2 y + y 3 ) + B 8I (x 3 + xy 2 ) + B 9I (x 2 + y 2 ) 2 = 0, (57d) 



where we have taken ( I = x + iy, the subscripts {IR} and {II}, 1 — 1,2, ...,9, represent 

the real and imaginary parts respectively. The expressions for Aj's and BiS are lengthy but 

can be obtained straightforwardly (by making use of (56) and the expressions (45a)) and so 

we do not present them here. Solving these overdetermined system of equations for x and 

y will give the suitable ratio I for which the shape restoring property of one of the 

V a 3 J 

solitons Si only arises in a three soliton collision process. 



T/ = l, T 2 ^l, Tf^l, j = 1,2. 



(56) 
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FIG. 7: Shape restoring property of soliton 1 (Si) during its collision with the other two solitons, 
soliton 2 (S2) and soliton 3 (S3), for the choice of the parameters k\ = 1 + i, k 2 = 1.5 — 0.5i, 
k 3 = 2 - i, af } = a< 2) = af = 1, = a « = 1.19, a f = =™ . 

Though we have not investigated the problem of existence of solutions of the equations 
(57), one can make a numerical search and identify suitable values of x and y to demonstrate 
the shape restoration property. For example in Fig. 7 with the parameters fixed at k\ = 
k 2 = 1.5 - 0.5*, h = 2 - i, a? = af ] = a< 2 > = 1, = ™, a« = 1.19, a?> = ™, 
we have demonstrated the shape restoration property. We find that while the amplitudes of 
the two of the solitons (S 2 and S3) change after interaction, the amplitude of the soliton Si 
remains unchanged during the interaction process. 

In the above analysis we have required the complete restoration property of soliton Si. 
However, it is also possible to require the intensity alone to be restored. In this case, the 
condition (56) can be modified as 

|T/| = 1, |T 2 |^1, |T?|^1, j = 1,2, (58) 

leading to a set of two complicated equations for x and y (which are too lengthy to be 
presented here). Solving them we can find x and y. Note that the quantities x and y 
correspond to the real and imaginary parts of the ratio of the parameters and af \ so 
that for every choice of x and y there exists a large set of and values for which 
shape restoration property holds good. 

One might also go a step further and demand that the phase shift $ 1 or the changes 
in the relative separation distances At i2 and At i3 vanish. These will give additional con- 
straints on the choice of parameters and These considerations require separate 
study and we have not pursued them here. It is obvious that such shape changing and 
shape restoring collision properties of the optical solitons in integrable CNLS equations, 
exhibiting a redistribution of intensity among the three solitons in the two modes, will have 
considerable technological applications both in optical communications including wavelength 
division multiplexing, optical switching devices, etc., and optical computation, for example 
in constructing logic gates [18, 19]. 
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G. Three-soliton solution of multicomponent CNLS equations and shape changing 
collisions 

The above analysis on the three soliton collision in 2-CNLS equations can be extended 
straightforwardly to the 3-soliton solution (14) of iV-CNLS equations, with arbitrary N, 
including N = 3. One can identify that shape changing collision occurs here also but 
with lot more possibilities for redistribution of intensities in contrast to the 2-CNLS case. 
The quantities characterizing the collision process here also are the intensity redistribution, 
amplitude dependent phase shifts and relative separation distances between the solitons, as 
explained in the 2-CNLS case. 

We also note that as the number of components increases from two to some arbitrary N 
(N > 2), the different possibilities for redistribution of intensity among them also increase 
in a manifold way. The corresponding transition matrix, measuring this redistribution, is 
found to be similar to Eqs. (45) with the redefinition of /%'s as given in Eq. (14b) along with 
the index j running from 1 to N instead of 1 to 2. The other factors, amplitude dependent 
phase shifts and change in relative separation distances, also bear the same form given by 
Eqs. (46) and (47), respectively, with this redefinition. 

As to the shape restoration property one has to again solve the equations 

T/ = l, Tf^O, Tf^O, j = 1,2, 3, ...,N. (59) 

Alternatively for intensity restoration the conditions are 

|T/| = 1, |T, 2 |^0, \Tf\^0, j = 1,2, 3, ...,7V. (60) 

Extending the above analysis, it is clear that, carrying out an asymptotic analysis of 4- 
soliton solution given in Appendix A, it is possible to restore the shape of two of the solitons 
at the maximum, which can be further generalized to arbitrary N soliton case, in which it 
is possible to restore the shape of N — 2 of the solitons. We have checked in this case also 
from the asymptotic analysis the soliton interaction is pair-wise and we conjecture that this 
should be true for the arbitrary A-soliton case as well. 
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VIII. MULTISOLITON SOLUTIONS AS LOGIC GATES 



The state vectors and LFTs introduced in Sec. VI and the shape changing pair-wise 
collision nature of bright solitons mentioned in Sec. VII can be profitably used to look at 
the multisoliton solutions of CNLS equations as various logic gates. We believe that such 
an approach provides an alternative point of view of shape changing soliton collisions to 
construct logic gates as discussed in Ref. [19]. The present point of view may have its 
own advantage as system initial conditions are chosen suitably to generate specific forms 
of multisolitons to represent logic gates may be much easier from a practical point of view, 
including replication, compared to constructing them through predetermined independent 
soliton collisions. In the following we will demonstrate this idea for the case of the 2-CNLS 
as an example. 



A. Three soliton solution and state restoration property 

The shape restoration of a particular soliton in arbitrary state associated with the three 
soliton solution has been discussed in Sec. VII F. Particularly this can be well appreciated 
with respect to binary logic states. For example, if we consider the soliton Si is in 1 state 
with the state value p\~ 2 — 1, it implies 

= (61) 

To obtain this we choose, ot^ = = 1. For simplicity we require S2 to be in state 

before interaction. From the asymptotic expressions (39), this can be achieved by choosing 

(i) 

the ratio %^ as 

h + kl 

a (2) 2k 2 + k{ - h ■ 1 ) 

Now in order to restore the state of Si after two collisions, we have to allow the outcome of 
Si resulting after the first collision, which may be called soliton S[, to interact with soliton 
S3 having a state inverse to the above state. This state for S3 can be identified from its 
asymptotic form before interaction given in Eq. (40). The resulting condition can be shown 
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to be 



where 



(i) 

a 3 n 



(63a) 



n = -(a? + 4 2) )a2 2) *^ + K 22 (h - k\ - 2k 3 )B + 1a?af*C 

-af{af* + af*)D + + a^\ 2 E, (63b) 
d = (a? + a { ^)a^*A - ^(^ + fc^)S - 2a^a?*C 

-(a^ + a^a^D, (63c) 

A = (h + k*)(k 3 + k*)(h + k*), (63d) 

B = (h + kDih + k^ih + k;), (63e) 

C = (h + kDih + kDih + k;), (63f) 

D = (k 2 + K)(k 3 + k*)(h + kl), (63g) 

£ = (fc 3 + k*)(h + kl)(k 3 + kl). (63h) 

In the above equation choosing the parameters satisfying conditions (61) and (62) one can 
fix ( ] suitably in order to restore the state of soliton Si. Thus the three soliton solution 



given by Eq. (10) having the specific choice of parameters specified by Eqs. (61)-(63) 
corresponds to the state restoration of soliton Si. 



B. Four-soliton solution and COPY gate 

Extending the above procedure, we can now consider the four-soliton solution given in 
Appendix A, and identify it as (i) a COPY gate or (ii) a ONE gate or (iii) a NOT gate 
studied in Ref. [19] for suitable choices of the arbitrary parameters. As an example, let us 
consider copying 1 state of Si to the output state of soliton S 4 . This requires the following 
steps. 

1. We consider the four soliton collision process in which the soliton Si collides with the 
soliton S 2 first and then with the soliton S3 and finally with the soliton S4. This 
sequence of collision follows from the condition ku > k 2 i > k 3 i > k 4 j. 

2. Consider for convinience Si to be in the state, the so-called actuator state [19]. This 
requires = 0, which can be obtained by choosing = and as arbitrary. 
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3. Assign 1 state to soliton £2 before interaction, for which we need 

(64) 



(dp k 2 — ki 



h + kf 

4. After its collision with S 2 as a result of shape changing collision the outcoming state 
of Si (say S[) will be altered. 

5. Now let us allow the third soliton in the four soliton solution to interact with S[ which 
changes the state S[ to S". 

6. Finally S4 is allowed to interact with S". From the asymptotic analysis, we identify 
the state of soliton S 4 after interaction as 



P\% = \y (65) 



(i) 



a. 



4 



We impose the condition on this state that this should be in the state of S2 before 
interaction. Thus the parameters and af^ of soliton S4 get fixed depending upon 
the input state of S 2 . 

7. The asymptotic analysis of the four soliton solution given in Appendix A results in 
following condition for S4 to be in one state after interaction. 

TtA\~ 1 



Tl At 



(66) 



where T-j 4 and are the transition elements of S4 in the modes qi and q 2 respectively. 
Here A^'k^R and A^k^R are the amplitudes of soliton S± before interaction in the two 
modes respectively. 

8. If we flip the input state of S 2 from 1 to state by suitably choosing the Pi^'s param- 
eters then the condition on soliton S4S output will become 

Pr4. = 0. (67) 

9. In the above two Eqs. (66) and (67) only free parameters are and af \ In principle, 
we can solve these two complex equations to obtain the free complex parameters 
and af 1 . Then for the given choice of parameters the state of the incoming soliton S 2 
can be copied on to the outgoing soliton S 4 . 
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Thus a four soliton collision process with the above premise is equivalent to a COPY gate. 
Similar procedure can be extended to other gates mentioned above as well. One can extend 
this idea further to identify a FANOUT gate from a 5-soliton solution. It appears that one 
can pursue the idea ultimately to identify the NAND gate itself as a multisoliton solution 
following the construction of Steiglitz in Ref. [19]. Fuller details will be reported elsewhere. 

IX. BRIGHT SOLITON SOLUTIONS AND PARTIALLY COHERENT SOLITONS 

As mentioned in the introduction, the recent observations by several authors [11, 12, 25] 
have shown that the Y-CNLS equations (1) can support Y-PCSs solutions. In general, 
these PCSs are said to be special cases of the so-called multisoliton complexes [2] which are 
nonlinear superposition of fundamental bright solitons. It has also been demonstrated that 
these PCSs are formed only if the number of components in Eq. (1) is equal to the number 
of solitons. Then it is quite natural to look for the 2-PCS, 3-PCS, 4-PCS, etc., as special 
cases of the two-soliton solution of the 2-CNLS, three-soliton solution of the 3-CNLS, four- 
soliton solution of the 4-CNLS equations etc., respectively, deduced in Sees. Ill and IV. In 
the following, we indeed show that the PCSs reported in Refs.fll, 12, 25] result as special 
cases, that is specific choices of some of the arbitrary complex parameters, from the bright 
soliton solutions of CNLS equations discussed in Sees. Ill and IV and thereby showing the 
origin of the various interesting properties of the PCS solutions. 

A. 2-PCS : A special case of the bright two-soliton solution of 2-CNLS equations 

Let us consider stationary limit of the two-soliton solution of the 2-CNLS equations 
(Manakov system) given by Eq. (8), that is, k n i = 0, for the special choice of the parameters, 
ocp = e mo , otp = —e V20 , = —otp = 0, where rj^s, j = 1,2, are now restricted as real 
constants. Then Eq. (8) becomes 




(68b) 



(68a) 
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FIG. 8: Typical 2-PCS forms for the Manakov system for z = see Eqs. (68), with k\ = 1.0 and 
k 2 = 2.0: (a) Symmetric case {t\ 2 = 0), (b) Asymmetric case (t\ 2 = 1). 



where 



D = 1 + n 



e m+vl e m+V2 

Ab 2 ' 



and 



+ 



\i 2 {k 



1R 



k 2R ) 2 e^ + ^ +m+ ^ 



16kl R k% R (k 1R + k 2R ) 2 



(68c) 



r}j = k jR (t + ik jR z) + Tfco, j = 1, 2. 



(68d) 



This stationary solution can be easily identified as the 2-PCS expression (13-15) given in 
Ref. [12] with the identification of x/s as t/s, j = 1, 2, 



t-t!=t+- + — — log 

h R 2k iR 

t-t 2 =t+- + — — log 

2fc 2 i? 



jj,(k 1R - k 2R ) 
Ak 2 1R (k 1R + k 2R ) 
n(ki R - k 2R ) 



(69a) 
(69b) 



_4:kj R (k 1R + k 2R )_ 

As the 2-PCS is a special case of the bright two-soliton solution of 2-CNLS equations, it is 
also characterized by a^'s (through rj^s) and fc^'s resulting in amplitude dependent phases 
and hence amplitude dependent relative separation distances. To be specific, in the PCSs the 
change in the relative separation distance play a predominant role in determining their shape 
as pointed out in [11,12]. These PCSs can be classified into two types as symmetric and 
asymmetric depending on the relative separation distances. Defining the relative separation 
distance i 12 = t 2 — t±, one can check that, for t 12 = 0, the PCS bears a symmetric form 
with respect to its propagation direction and is known as symmetric PCS [11]. It takes an 
asymmetric form for t\ 2 ^ and is known as asymmetric PCS [11]. From Eqs. (69), the 
relative separation distances for the stationary 2-PCS can be obtained as 



t 12 = t 2 -t 1 = - — + — — log 

KlR K 2R z.Kl R 



/i(k 



1R 



k 



2 R ) 



4=k 2 R (k 1R + k 2R ) 



2k 



2 R 



-log 



fj,(k 



1R 



k 



2R) 



4k 2 R (k 1R + k 2R ) 



(70) 



Typical forms of symmetric and asymmetric stationary 2-PCS are shown in Fig. 8, which 
in similar to the ones in Ref. [12]. 



39 



B. 3-PCS: A special case of the bright three-soliton solution of 3-CNLS equations 



Since it has been observed that the PCS solutions exist when the number of components 
is equal to the number of solitons propagating in the system, we consider next the 3-soliton 
solution of the 3-CNLS equations in order to show that the 3-PCS as a special case of the 
3-soliton solution. Thus considering the stationary limit k n i = 0, n = 1, 2, 3, of the 3-soliton 
solution of 3-CNLS equations given by Eq. (14) with N = 3, and making the following 
parametric choice 



e mo ,a { 2 ] 



-e" 2 °,c4 3) 



(3) 



4" 



(3) 
«2 



a 



(i) 



a 



(2) 



0, (71) 



where r/jo's, j = 1, 2, 3, are restricted to real parameters, we obtain 



<?2 = 



Q3 



Here, 



+ 



v{kiR-k 2R )e^2+v* 2 n{k 1R -k 3R )e^ +r * + ^ 



+ 



4k 2 R {k 1R + k 2R ) *kl R {k m + k 3R ) 

| fi 2 (k 2R - k 1R )(k 3R - k 1R )(k 3R - fc 2R )V3+^ 2 +^+rn 
l6k 2R k m( k 2R + kiR>(h R + k 1R )(k 3R + k 2R ) 2 
oV2 Kkm - k 2R )e^ + ^ fx(k 3R - k 2R )e«* +T 'Z + ** 
^k 2 R (k lR + k 2R ) 4:k 2 R (k 3R + k 2R ) 

, H 2 (k 2R - k 1R )(k 3R - k 2R ){k 3R - k 1R ) 2 e^ + ^ + ^ 



oV3 



+ 



\Qk 2 lR kl R {k 2R + ki R ){k 3R + k 2R ){k 3R + k 1R ) 
Kk 3R -k 1R )e^+V3 n{k 3R -k 2R )e^ + ™ 



4k 2 R {k 1R + k 3R ) 



+ 



^ k m( k 3R + k 2R ) 



V 2 (k 3R - k 1R ){k 3R - k 2R )(k 2R - fc lfi ) 

16k 2 R k 2 R (k 3R + k 1R )(k 3R + k 2R )(k 2R + k lR ) 2 



D 1 = l+/i 



4*^ 



+ 



4%R 



+ 



M 2 3R 



+ 



H 2 {k 



/D 



/A 



IR 



^2R 



16k 2 R k 2 R (k 1R + k 2R ) 



. f^iR ~ k 3R ) 2 e^l+m+nt _ k2R)2eV2+v * +m+v * 3 



+ 



1Qk lR k 3R( k lR + k 3R? 1Gk 2R k 3R( k 2R + k 3R? 

H 3 (k 2R - k 1R ) 2 (k 3R - k 1R ) 2 (k 3R - k 2R ) 2 e^ + ^ + ^ 

^ k lR k 2R k 3R( k lR + k 2R) 2 ( k lR + k 3Ri) 2 {. k 2R + k 3R) 2 



(72a) 



(72b) 



(72c) 



(72d) 



The above solution can be easily rewritten as Eq. (16-18) for the 3-PCS case given in Ref. 
[12]. As in the case of 2-PCS, here also we identify %'s given in Ref. [12] as t/s, j = 1, 2, 3, 
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which are defined below: 



t-u 



t-u 



km ZkiR 



t + 



V20 



1 



k 



2R 



2k 



2R 



-log 



K 3R ^K 3 r 



fj,{k 2 R - k 1R )(k 3R - km) 
Akl R (k m + k 2R ){k 1R + k 3R ) 

n(k 2R - k 1R )(k 3R - k 2 R) 
Akl R (k 1R + k 2R )(k 2R + k 3R ) 

H\km — km)(k 3 R — k 2R ) 



(73a) 
(73b) 
(73c) 



_4kl R (k 1R + k 3R )(k 2R + k 3R ) 
These 3-PCSs can also be classified as symmetric and asymmetric as in the case of 2-PCSs. 
The stationary 3-PCS is symmetric when t 12 — t 13 — and asymmetric otherwise. In Fig. 9 
we have shown the symmetric and asymmetric 3-PCS solutions. 

FIG. 9: Typical 3-PCS forms for the integrable 3-CNLS system for z = with fcj = 1.0, k 2 = 0.5 and 
k 3 = 0.2, see Eqs. (72): (a) Symmetric case (ti 2 = t\ 3 = 0), (b) Asymmetric case (t\ 2 = 1, t\ 3 = 2). 



C. 4-PCS: A special case of the four-soliton solution of 4-CNLS equations 



In a similar fashion as in the above two cases, the 4-soliton solution of the 4-CNLS 
equations given in Appendix A, with = 4 can also be shown to reduce to 4-PCS given 



by Eq. (19-23) in Ref. [12] by choosing, k nI = 0, a\ 



(i) 



(2) 



— pV20 



a 



(3) 



a 



(4) 



e mo , at*p = 0, i,j = 1, 2, 3, 4, i ^ j. Since it is straightforward but lengthy to write 



down the form, we desist from presenting the solution here. Here the t,'s, j = 1,2,3,4, are 
defined as 



to = 
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fi(k 2R - k 1R )(k 3 R - k 2R )(k m - k 2R ) 
^k\ R {k 2R + k 1R )(k 3 R + k 2R )(k 4R + k 2R ) 

n{k 3 R - km)(k 3R - k 2R )(k AR - k 3R ) 
Akl R (k 3R + km)(k 3R + k 2R )(k AR + k 3R ) _ 

H{k m - k 1R )(km ~ k 2R )(k 4R - k 3R ) 



(74a) 
(74b) 
(74c) 
(74d) 



_4kl R (k 4R + km){k 4 R + k 2R )(k 4R + k 3R ) _ 
Here also the symmetric PCS results for = 0, j > i, i,j = 1, 2, 3, 4, and asymmetric PCS 
for Uj ^ 0, j > i. 

Extending this idea to arbitrary N, it is clear that the A^-PCS is a special case of the 
Af-soliton solution of A^-CNLS equations (1). It has been noticed in Refs. [11, 12] that these 
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PCSs are of variable shape. The reason for the variable shape can be traced naturally to 
the nontrivial dependence of phases on the complex parameters a^'s as shown above. Thus 
it is clear that any change in the amplitude will affect the phase part of the solitons and 
vice- versa. Since we have explicitly shown that iV-PCSs are special cases of bright N-soliton 
solutions of iV-CNLS equations, they possess variable shape as a consequence of the shape 
dependence on the af^ parameters. 

D. Propagation of partially coherent solitons and their collision properties 

The intriguing collision properties of the partially coherent solitons reported in Refs. [11, 
12] can be well understood by writing down the expression for PCSs with nonvanishing Aw's, 
that is nonstationary special cases of multicomponent higher order bright soliton solutions 
discussed in Sees. Ill and IV. For the nonstationary PCSs we can choose as a special case the 
complex parameters a-^'s (i ^ j) to be functions of fc„/s such that they vanish as k n i = 0. 
As we make these k n j ^ 0, then a^'s (i ^ j) also vary, thereby making the collision scenario 
interesting. We can consider both the cases of equal and unequal velocities, which exhibit 
similar behaviours. 

As a first example, we consider the propagation of the 2-PCS comprising two solitons 
with equal velocities (ku = k 2 i) in PR media. Its propagation can be studied by choosing 
(for illustrative purpose) = ku and = (0.25 + 1.02i)fc2j as functions of velocities 
= 1,2) such that they vanish when kji = 0, j = 1,2. This is shown in Fig. 10 for 
the parameters ol^ = 2.0 + i, = 1, ki — 1.0 + i, and k 2 = 2.0 + i. For unequal 

velocity case (ku ^ fc 2 /), the PCS collision is shown in Fig. 11 for the parametric choice 
a? = 1.0, af ] = k u , = -{^^ L )k 2 i, ^ = -2.0, k x = 1.0 + i, and k 2 = 2.0 -i. This 
can also be viewed as collision of two 1-PCS which are spread among the two components, 
which are travelling with equal but opposite velocities. 

FIG. 10: Intensity profiles showing the collision scenario of two 1-PCSs, with equal velocities, at 
(a) z = — 5 and (b) z = 5, given by special choice of parameters (as given in text) in the 2-soliton 
solution of the Manakov system. 

Now let us consider the collision of 2-PCS and 1-PCS in PR media. This is equivalent 
to the three soliton collision in the 3-CNLS system with specific choice of parameters. We 

42 



FIG. 11: Intensity profiles showing the collision scenario of two 1-PCSs, moving with equal but 
opposite velocities, at (a) z = —5 and (b) z = 5, given by special choice of parameters (as given in 
text) in the 2-soliton solution of the Manakov system. 

consider the case, in which the complex parameters Oi[\ af^ , , a 2 3 \ ol%\ a^ 1 are 
nonvanishing and as functions of fc n j's, n = 1,2,3. Then the resulting asymptotic forms of 
the 3-PCS propagation is shown in Fig. 12 for the parametric choice cvf^ = 1.0, 
a f) = k u , 0$ = -0.5k 2I , af ] = 0.25, af = 0.02A; 2/ , «f' = ~(^w^)hi, af = 2 A*/, 

(3) 

a 3 = —2, ki = 1.0 + i, k 2 — 1.5 — i, and fc 3 = 2.0 — i . In the above figures it can be 

FIG. 12: Intensity profiles showing the collision scenario of 2-PCS with 1-PCS at (a) z = — 4 and 
(b) z = 4 given by special choice of parameters (as given in text) in the 3-soliton solution of Eq.(l) 
with N = 3. 

verified that the total intensity of the individual solitons comprising the PCS is conserved. 

The above analysis on PCS propagation clearly shows that, there will be variation in 
the shape of the PCS during its collision with other PCSs. The explanation for this result 
follows from the shape changing (intensity redistribution) nature of fundamental bright 
soliton collision of the integrable CNLS equations, explained in Section V. Further, we have 
also observed that the collision of two PCSs each comprising m and n soliton complexes, 
respectively, such that m+n = N studied in Refs. [11, 12, 25], is equivalent to the interaction 
of N fundamental bright solitons (for suitable specific choice of parameters) represented by 
the special case of iV-soliton solution of the iV-CNLS system. It should also be noted that 
in the collision process the total intensity of individual solitons comprising the iV-PCS is 
conserved. This is due to the complete integrable nature of the iV-CNLS equations (1). 

E. Multisoliton Complexes 

In the above we have considered the CNLS equations with number of components (say 
p) is equal to the number of fundamental solitons (say q). This is only a special case of 
the multisoliton complexes and has been much discusssed recently. However, the results are 
scarce for the case p ^ q, except for the work of Sukhorukov and Akhmediev [25], where 
the incoherent soliton collision is demonstrated numerically. To elucidate the understanding 
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we present a form of the three soliton complex in which three solitons are spread among 
the two components, by suitably choosing the parameters in the explicit expression, Eq. 
(10). This has been shown in Fig. 13 with the parameters chosen as = af 1 = 1.0, 
af = 0.5, c4 2) = 0.25, = af = -2, k x = 1.0 + i, and k 2 = 1.5 + i. From 

FIG. 13: Intensity profiles of a multisoliton complex comprising three solitons spread up in two 
components propagating in photorefractive media: A special case of the three soliton solution (10) 
of the integrable 2-CNLS system for the parameters chosen as in the text, (a) at z = — 5 and (b) 
at z = 5. 

the figure and the analysis of the soliton interaction it is clear that the shape of these 
complexes strongly dependent on the ct^'s along with k/s which determine how the solitons 
are spread up among the components. For the same case there exists various forms of 
multisoliton complexes depending on the spreading up of solitons in the two components. 
As a consequence of this multisoliton complexes will possess a rich variety of structures in 
comparison with the PCSs. 

X. CONCLUSION 

We conclude this paper by stating that the collision processes of solitons in coupled 
nonlinear Schrodinger equations lead to very many exciting novel properties and potential 
applications. The novel properties include shape changing intensity redistributions, ampli- 
tude dependent phase shifts and relative separation distances, within the pair-wise collision 
mechanism of solitons. Interestingly, it is identified that the intensity redistribution char- 
acterising the shape changing collision process in A^-CNLS equations can be written as a 
generalized linear fraction transformation. This will giver further impetus in constructing 
multistate logic, multi-input logic gates, memory storage devices and so on by using soliton 
interactions. The implication of these properties require further deep investigations. Fur- 
ther, viewing the recently much discussed objects multisoliton complexes, partially coherent 
solitons as special cases of the bright soliton solution enhances the understanding of their 
various properties. We expect the interaction study presented here will shine more light on 
spatial soliton propagation in (1 + 1)D photorefractive planar waveguides. 
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APPENDIX A: FOUR-SOLITON SOLUTION 

In this appendix for completeness, we present the form of the four-soliton solution of the 
2-CNLS equations by generalizing the two- and three- soliton solutions of it, which can be 
obtained by terminating the power series as 



and solving the resulting set of linear partial differential equations recursively. It can be 
written as 



9 U) = xg? + x 3 g? + x 5 g? + x 7 g?, 

f = 1 + X 2 f2 + X 4 U + X 6 fe + x7s, 



(Ala) 
(Alb) 
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In the above 
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FIGURE CAPTIONS 



Fig. 1: Two distinct possibilities of the shape changing two soliton collision in the 
integrable 2-CNLS system. The parameters are chosen as (a) k\ = 1 + i, k 2 = 2 — i, 
of = af = af = 1, af = 22±p; (b) k x = 1 + i, k 2 = 2 - i, of = 0.02 + O.li, 
af = af = af = 1. 

Fig. 2: Intensity profiles of the three modes of the two-soliton solution in a waveguide 
described by the 3-CNLS (Eq. (1) with N = 3) showing different dramatic scenarios of the 
shape changing collision for various choices of parameters. 

Fig. 3 : Plot of the magnitude of phase shift as a function of the parameter af , when 
it is real (for illustrative purpose), see Eqs. (29-31). The other parameters are chosen as 
fci = 1 + %, k 2 = 2 - %, af = af = 1 and q?' = 

Fig. 4: A schematic three soliton collision process (for the choice km, k 2 R, k 3 R > 0, 
ku > k 2I > k 3I ). The effect of phase shifts are not included in the figure. 

Fig. 5: Intensity profiles \qi\ 2 and \q 2 1 2 of the two modes of the three-soliton solution of the 
2-CNLS equations, representing elastic collision, with the parameters chosen as k± — 1 + i, 
k 2 = 1.5 — 0.5i, k 3 = 2 — i, ol^ = a 2 ^ = ajp = = a 2 ^ = = 1. 

Fig. 6: Intensity profiles |gi| 2 and |g 2 | 2 of the two modes of the three-soliton solution of 
the 2-CNLS equations, representing the shape changing (intensity redistribution) collision 
process for the choice of the parameters, k\ — 1 + i, k 2 — 1.5 — 0.5i, k$ = 2 — i, = 39 ~ 9 80i , 

a (l) = 39^^ a (l) = Q g + Q ^ a (2) = Q a (2) = a (2) = L 

Fig. 7: Shape restoring property of soliton 1 (S'i) during its collision with the other 
two solitons, soliton 2 (S 2 ) and soliton 3 (S3), for the choice of the parameters k± — 1 + i, 
k 2 = 1.5 - 0.52, k 3 = 2 - 1, a? = af ] = a { 2 2) = 1, = ™, af = 1.19, af = ^ . 

Fig. 8: Typical 2-PCS forms for the Manakov system for z = see Eqs. (68), with 
ki — 1.0 and k 2 = 2.0: (a) Symmetric case (t 12 = 0), (b) Asymmetric case (t 12 = 1). 

Fig. 9: Typical 3-PCS forms for the integrable 3-CNLS system for z = with k\ = 1.0, 
k 2 = 0.5 and k 3 = 0.2, see Eqs. (72): (a) Symmetric case (t± 2 = tis = 0), (b) Asymmetric 
case (t 12 = 1, t 13 = 2). 

Fig. 10: Intensity profiles showing the collision scenario of two 1-PCSs, with equal 
velocities, at (a) z = —5 and (b) z — 5, given by special choice of parameters (as given in 
text) in the 2-soliton solution of the Manakov system. 
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Fig. 11: Intensity profiles showing the collision scenario of two 1-PCSs, moving with equal 
but opposite velocities, at (a) z = — 5 and (b) z — 5, given by special choice of parameters 
(as given in text) in the 2-soliton solution of the Manakov system. 

Fig. 12: Intensity profiles showing the collision scenario of 2-PCS with 1-PCS at (a) 
z = —4 and (b) z — 4 given by special choice of parameters (as given in text) in the 
3-soliton solution of Eq.(l) with N — 3. 

Fig. 13: Intensity profiles of a multisoliton complex comprising three solitons spread up 
in two components propagating in photorefractive media: A special case of the three soliton 
solution (10) of the integrable 2-CNLS system for the parameters chosen as in the text, (a) 
at z = — 5 and (b) at z — 5. 
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